Transport on percolation clusters with power-law distributed bond strengths: when 

do blobs matter? 
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The simplest transport problem, namely maxflow, is investigated on critical percolation clus- 
ters in two and three dimensions, using a combination of extremal statistics arguments and exact 
numerical computations, for power-law distributed bond strengths of the type P(cr) ~ o~ a . As- 
suming that only cutting bonds determine the flow, the maxflow critical exponent v is found to 
be v(a) = (d — l)v + 1/(1 — a). This prediction is confirmed with excellent accuracy using large- 
scale numerical simulation in two and three dimensions. However, in the region of anomalous bond 
capacity distributions (0 < a < 1) we demonstrate that, due to cluster-structure fluctuations, it 
is not the cutting bonds but the blobs that set the transport properties of the backbone. This 
"blob-dominance" avoids a cross-over to a regime where structural details, the distribution of the 
number of red or cutting bonds, would set the scaling. The restored scaling exponents however still 
follow the simplistic red bond estimate. This is argued to be due to the existence of a hierarchy 
of so-called minimum cut-configurations, for which cutting bonds form the lowest level, and whose 
transport properties scale all in the same way. We point out the relevance of our findings to other 
scalar transport problems (i.e. conductivity). 

PACS numbers: 64.60.Ak,74.25.Fy,73.50.-h 



I. INTRODUCTION 



The transport properties of percolation clusters have 
been a subject of interest for many years g|. A 
natural problem to study is e.g. conductivity, and 
one often complicates it further by using random bond 
"strengths" a with a power-law tail of the form a~ a 
{§ § § § 0, |, In the first place, because this 

allows to represent continuum percolation |5|, Q and thus 
get closer to some actual physical realizations of percola- 
tion. A second reason why these systems are interesting 
is the equivalence |l2j between transport on strongly dis- 
ordered systems and percolative transport ]35| . 
Transport critical exponents on these systems are found 
to depend on a, which means that strict universality is 
lost. The original observation that transport exponents 
become nonuniversal is due to Kogut and Straley [3 
who used mean- field type arguments. Later Straley [3 
with the help of the nodes-links-blobs [HJ, [l5| |l_ 
picture of the backbone, concluded that the conductivity 
exponent t, such that £ ~ (p — p c ) , is the maximum 
of the universal exponent to — (d — 2)v + £ and the in- 
dependent exponent t(a) = (d— 2)^+1/(1 — a). Here v is 
the correlation length exponent and ( measures the con- 
tribution of blobs to the resistance between two points 
on the backbone, for the case of constant conductances. 
For the conductivity problem thus there is a crossover 
from the universal exponent to for a < a c to t(a) in the 
"anomalous regime" a > a c . Although not without some 
controversy initially M, H , ^, [l7j] , this result is by now well 
established §, g, 0, yT 



It is somehow surprising that i{a) can be analytically 
calculated in the "anomalous regime" , given that the uni- 
versal exponent to, which applies to the arguably simpler 
case of constant conductance, has not been analytically 
derived up to now. The difficulty in deriving to resides 
in that £ is determined by the blobs, and thus one would 
require detailed information |lS| ] about the internal struc- 
ture of the blobs. 

On the contrary, it has been argued by several au- 
thors H [ll], that t(a) in the anomalous regime is 
determined by the cutting-bonds alone. Since these form 
linear chains of typically l}l v bonds at p c jL5|, the re- 
sulting conductivity exponent is easily derived. The ar- 
gument to support the belief that blobs are irrelevant in 
the anomalous regime seems to be roughly the following: 
an exceedingly small conductivity falling on a blob has 
little effect on the overall conductance, because there are 
many alternative parallel paths. On the other hand, if 
this small conductivity is located on a cutting bond it 
will certainly dominate the system conductance. While 
this argument is true in principle, this reasoning misses 
the fact that the number of cutting bonds is itself a fluc- 
tuating quantity. The issue of blob-irrelevance has been 
considered by Machta and collaborators Q using a hi- 
erarchical model for the backbone, to reach similar con- 
clusions. However, as noted by the authors, their model 
does not include structural fluctuations. We will demon- 
strate in this work that it is in fact the blobs and not 
the cutting bonds, that determine the critical transport 
properties, even in the anomalous regime. This as we 
will see, is due to structural fluctuations. However, the 
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resulting transport properties turn out to be the same as 
those given by the most simplistic red bond estimate. 

A related critical transport problem, which is relevant 
for disordered superconductors, is that of determining 
the critical current density J c = l c L^^ d ^ that a perco- 
lation network can sustain, and which above p c behaves 
as J c ~ (p - Pc ) v |J, || §, i||p. This prob- 
lem has a simple geometrical interpretation. Finding the 
maximum flow of current, or maxflow, is equivalent to 
finding the surface across the system, on which the sum 
of critical currents of the bonds is maximized. As we 
will draw advantage of this analogy later on, we note 
that this surface is called a mincut in computer science 
language (2§[§,[27). 

In this paper our goal is to present the first comprehen- 
sive study of the maxflow problem on percolation clus- 
ters. This is motivated by the following observations. 
First, this is the simplest transport problem that one can 
think about, and has not been as such discussed much in 
the literature. Second, we are able to use to our advan- 
tage recent developments [p8| on combinatorial optimiza- 
tion algorithms, in the context of disordered systems. 
Here one can use a three-step approach, in which first 
a critical spanning cluster is set, its backbone is pruned 
out, and finally that is used for the max-flow/min-cut 
problem. Each stage is solved with one of the powerful 
graph optimization algorithms for the particular prob- 
lem, as discussed later. 

In the simplest version of the maxflow problem, all 
present bonds have the same critical current, or capacity 
i c and absent bonds have i c = 0. At criticality, a typical 
percolating cluster is a linear chain of cutting bonds and 
thus I c = i c - From this observation plus the usual scaling 
relation J C (L) ~ L~ v l v , one concludes that v = (d — 
l)v. This result is consistent with experiments [[l9[ |2^1 
and numerical simulation j^], In a more realistic 

model, each present bond has a random capacity i c with 
power-law distribution P(i c ) ~ i^ a - This is for example 
the case for continuum percolation models (20[ [H], p3| . 
A simple extension of the "typical cutting bond string" 
argument gives v(a) = (d — \)v + 1/(1 — a) as we show 
later in Section y. 

In the following we will find it useful to compare 
the conductivity and critical current problems to each 
other. This comparison is done by interpreting the ran- 
dom bond variables i c alternatively as bond conductances 
a or as bond capacities i c . Consider for example two 
bonds with o\ and a 2 connected in parallel. The re- 
sulting conductance a par = o~\ + a 2 is then the same as 
the maximum current I ma x that can flow if o~j are ca- 
pacities. If these bonds are instead connected in series, 
then aseries = (cf 1 + v^y 1 and I max = min(ai,a 2 ) 
are no longer equal. However the series conductance can 
be written as p9| a ser i es — min{a\, 02XI + P) _1 j with 
/3 = min(cri, G2)Jrnax{(j\, 02)- In the limit of strong dis- 
order (a — > 1), is typically negligible. We conclude 
that, in this limit, also in the series case the conductance 
equals exactly the maximum current obtained by inter- 



preting <7j as capacities i c . Therefore in the a — > 1 limit, 
the resistive current problem and the superconducting 
current problem (maxflow) are equivalent, at least for all 
structures that can be solved by a combination of series 
and par allel b ond reductions ]36[ . Moreover as shown in 
Section [II B , we find that the equivalence noticed above 
is valid not only in the a — > 1 limit but for a range of a 
values, for strings of bonds in series. 

In deriving a-dependent exponents, both for v(a) and 
i(a), the assumption is made that the backbone always 
contains L x l v cutting bonds. While this is true typi- 
cally, the number of cutting bonds is in fact a fluctuating 
variable whose distribution may extend down to zero in 
the form of a power law (see later). The existence of 
such fluctuations has been noted by some works previ- 
ously , but their role in transport properties has 
not been considered. These number fluctuations, we will 
show in Section H, do modify the transport exponent 
that results from a string of cutting bonds. Then by an- 
alyzing the conceptually and numerically simple maxflow 
problem, we will be able to show that in fact blobs can- 
not be neglected. The net outcome, which we justify 
by a heuristic hierarchical picture, is that although the 
simplest cutting bond scaling (without fluctuations) is 
restored, it is in fact the blobs that set this scaling be- 
havior. 

The structure of the rest of the paper is as follows. 
Section || presents the analytical discussion, based on a 
"fluctuating number of cutting-bonds" picture. In Sec- 
tion |H we go through one-by-one the numerical methods 
employed, the findings about structural fluctuations, and 
some further numerical analysis of the extremal statistics 
aspects. Section [lV| contains the results concerning the 
maxflow problem, and some details of interest that can 
be determined from analyzing large statistics. Section |v| 
finishes the paper with a discussion. 



II. CRITICAL CURRENT DENSITY 

We consider diluted lattices where the maximum su- 
percurrent i c that a present bond can sustain is a random 
variable distributed between and I according to 



P(i c ) = (1 - a)c~ a , 



(1) 



with a < 1. 

Let I c be the maximum supercurrent (or maxflow) 
that the whole system, given a set of values {i c }, 
can sustain. The average current density J c is then 
J c =< I c > I L^ -1 ), and goes to zero at p c as 



Jc~(p-Pc)" 



(2) 



Right at Pc, and for a system of finite linear size L, usual 
finite-size scaling arguments (TR Q imply that 



J c (Pc,L)~L- v ' v , 



(3) 
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where v is the percolation correlation length exponent. 
The nodes-links-blobs picture of the percolation clus- 
ter jl3|, [l4|, [ll| [l6| tells us that, right at p c , there is 
typically a single connected path through the sample. 
This path is a sequence of multiply connected regions 
(blobs) connected by strings of singly connected bonds, 
also called cutting bonds. The average number of cutting 
bonds is of order L 1 /" at p c fUsf . 

We now start by considering the maximum flow /* al- 
lowed by a string of n cutting bonds, and which obviously 
equals the least capacity among the n bonds. The typi- 
cal least value /* among a collection of n » 1 random 
numbers i c with probability P(i c ) satisfies 



/n 



Thus 



f 

J r 



P{i)di = 1/n 



= „-V(l-a) 



(4) 



(5) 



On a system of linear size L at p Cl the average number of 
cutting bonds is L 1 ^ mM. In replacing this one obtains 
/* - L-iMi-") _ LF^ 1 J C ~ L(d-i)-v(a)/^ and thus 

from Eq. fit), 



v(a) = (d-l)i/+l/(l-a) 



(6) 



as advanced in the introduction. 

This typical-n argument however neglects the fact that 
n is a fluctuating number. Since V(ji) actually has a 
power law tail extending down to n — ]3l| , this neglect 
turns out not to be correct for quantities that depend on 
1/n as Eq. (|). 

We now present a more careful treatment, which takes 
into account the fluctuations in n. It is known pl| ] 
that Vl(ii) = (n^) _1 "P(rt/n^), where V(h) is a size- 
independent function, and n* L ~ L 1 /" |l5|j . Since for 
the purpose of our discussion all that matters is the 
behavior of V(n) as n — ► 0, we take for simplicity 
V[h) = (1 + a)n a , for < n < 1. Thus, 



7>(n) = (l + o)(n£) 



-(l+a) n a 



(7) 



for 1 < n < n.^. We will for the moment assume that n 
cannot be zero. 

Let now / be the minimum among n numbers x dis- 
tributed with probability P[x). The distribution m n (f) 
of / is determined as 

m n (f) = nP(/)il- [ P(x)dx\ 



r>l f \ ~ n I P(x)dx 

iP(f)e J o . 



(8) 



Because of the strong exponential suppression that oc- 
curs for / larger than /* defined by Eq. (|j), we can ap- 
proximate m n (f) by 




if < / < /* 
if / > /* 



(9) 



Now allowing for the fact that n fluctuates, the probabil- 
ity distribution function (PDF) of the maxflow / through 
a string of cutting bonds is 



m{f) 



driP[ri)m n (f) 



= (i + «)K) 



-(i- 



dnn a m n (f), (10) 



for < / < 1. From Eq. @ and Eq. @ we conclude 
that, for a given value of /, the only nonzero contribu- 
tions in Eq. ( pL0| ) come from n values which are smaller 
than r}(f) = f^ 1 "^ 



Thus 



P(f)(\+a) r min «Mf)) 
m ff) = n />\ L J a) / n a+1 dn (11) 

Defining f typ = (n2)" 1/(1 " a) , « = (a + 1)(1 - a) and 
A = k/(cl + 2), this last expression can be written as 



x &yp (fL) 



if < / < /, 



' yp 



if ftyp < / < 1- 



(12) 

This gives the PDF for the maxflow / through a string 
of cutting bonds on a system of size L, allowing for fluc- 
tuations in the number n of bonds on the string. The 
strength of the fluctuations of 1/n is characterized by 
the exponent k, which in turn depends on a. If a — > oo 
(non-fluctuating limit) m(f) is nonzero only for / < ftyp- 
Thus < / >~ ftyp and Eq. (||) is recovered in this case. 
However it is known that a rj 0.22 in two dimensions |jl|| . 

For general a and a, m(f) has a power-law tail with 
exponent (1 + k) for f » ftyp- The importance of this 
power-law tail is evidenced by considering the average 
flow 



</> 



A 



A 



2 - a 



ftyp + {ftypT (13) 

_L " ri 

When k > 1, < / >~ ftyp - (nl)- 1 ^ 1 -^ - L -1 M 1-a 
and Eq. (||) is recovered. However if a > a /(a + 1) (re < 
1), the power-law tail dominates the average. In this case 
< / > ftyp » ftyp- Therefore < / >~ L-( a + 1 )/ y , an d 
Eq. @ implies that in this case, 



(d-l)u + a+l 



(14) 



The meaning of this is clear. If a is large, typical cases 
with 0(n* L ) cutting bonds will only allow an exceedingly 
small flux /. The average flow < / > however will be 
dominated by the very rare cases in which n is small and 
for which / ~ 0(1) >> ftyp- So finally we conclude that, 
if we idealize the backbone at p c as a string of n cutting 
bonds, and if V(n) behaves for small n as n a , one has 
that 

v(a) = { (15) 
a + 1 + v(d- 1) if a > ^ 
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III. NUMERICAL RESULTS 
A. Algorithms 

In this section we test our analytical derivation of 
v(a) of Section [il] in two and three dimensions on large 
systems, with the help of powerful combinatorial algo- 
rithms |^8| . Percolation backbones are first generated by 
means of a matching algorithm |p8| , [33[ , for square and 
cubic lattices. We do this by randomly adding bonds one 
at a time until a percolation path is first present. At this 
point the matching algorithm identifies the conducting 
backbone, exactly at the percolation point for each sam- 
ple. Alternatively one could fix the density of present 
bonds to a value close to the infinite system critical den- 
sity p c , and then identify the percolating backbone with 
the same algorithm. However, our procedure has the ad- 
vantage that no separate estimate is necessary for p c . 

For each percolating backbone, capacities are drawn 
from the given distribution, and the maxflow is calcu- 
lated by means of a flow augmentation algorithm (see 
p8[ for a review of the max- flow problem). The effi- 
ciency of the maxflow algorithm is highly increased when 
working on the backbone only, so we are able to analyze 
thousands of samples for each value of a. In this way we 
estimate numerically the average flow at p c for several 
linear sizes L, and from its scaling properties v(a) is de- 
rived. The largest sample sizes studied were L = 4000 in 
two dimensions and L = 120 in three dimensions. These 
are mostly set by the CPU usage of the combination of 
the matching and flow algorithms, which in turn is dom- 
inated for L large by the scaling of the matching part. 
The maxflow code is actually sub-linear in n = L 2 in 
CPU time, since the mass of the backbone scales with 
its fractal dimension. Notice that once the backbone of 
a sample has been established, it can be used for sev- 
eral consequent maxflow determinations for different a 
to save CPU-time. In the Appendix we present an idea 
for an optimal algorithm for this problem. 



B. Results 

Results are shown in Fig. [l]. Our numerical simulation 
results confirm Eq. (^) nicely. However the saturation 
of v(a) predicted by Eq. ( plf ) for a > a /(a + 1) does 
not occur. Notice that a is not a universal exponent but 
depends on the ensemble. For example if the ensemble 
is determined by fixingp = p c , numerical measurements 
and RG calculations |3l[] give a ~ 0.22. Additionally, 
roughly 20% of the connected samples have zero cutting 
bonds §(], |1J, that is, V(n) ~ 0.20S(n) + cn°- 22 for small 
n. 

However other ensembles can be considered. Consider 
for example the percolation cluster defined by Ambe- 
gaokar et al construction, in which conductances are laid 
down on the lattice in order of increasing conductivity un- 
til a percolating path is created |L2| . At least the last con- 
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FIG. 1: Maxflow exponents measured for 2d (empty squares) 
and 3d (filled squares) percolation clusters at p c . Dotted lines 
indicate the theoretical result v(a)/v — (d — 1) = 1/^(1 — ct), 
with 1/v = 0.75 (2D) and 1.13 (3D). The maximum linear 
size simulated was L — 4000 in 2d and L = 120 in 3d. Error 
bars are smaller than symbol sizes. 



ductance to be laid down is a cutting bond, so 'P(O) = 0. 
Experimentally this situation is realized when supercon- 
ductive samples grow percolatively by deposition |2j| . In 
this case the point at which the supercurrent is nonzero 
for the first time is defined by the first appearance of a 
connected path, not by a fixed density of occupied bonds. 

As we add bonds one at a time, our numerical simula- 
tions correspond to this case rather than to fixing p = p c . 
Our measurements of the distribution of the number of 
cutting bonds indicate (Fig. ||) that V(n) ~ n a for small 
n, with a w 1.25 in two and three dimensions. 

Different ensembles give rise to different distributions 
of the number of cutting bonds, and specifically to dif- 
ferent values of a so, if Eq. ([IB]) were to hold for perco- 
lation clusters, the resulting transport exponent would 
be ensemble- dependent. However, our maxflow measure- 
ments on percolation clusters are consistent with Eq. (|^) 
for all a, without any sign of saturation. 

In view of the failure of percolation clusters to show the 
predicted exponent saturation, we first confirmed the va- 
lidity of Eq. (15) for strings of cutting-bonds. We did 
so by numerically studying strings of bonds whose num- 
ber n is distributed according to (Q), and whose con- 
ductances (or capacities, for the maxflow problem) are 
distributed according to (|l]). The maximum flow is sim- 
ply the least critical current i c . Alternatively, bond ca- 
pacities i c may be interpreted as conductances, in which 
case the resulting conductance for the whole string is 



simply a = 1/X^=i 1 Ac (:?')■ We find for these strings 
of cutting bonds (Fig. ||) that Eq. ( |l5| ) is satisfied very 
accurately. Fig. || also shows that the conductivity and 
maxflow exponents are the same for a > 0, indicating 
that the conductance is dominated by the least i c value 
in that regime. We conclude that Eq. (fL5|) is exact for 
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FIG. 2: Distribution of the number n of cutting 
percolation clusters in two (top: L — 32 (empty 
128 (asterisks) and 1024 (full squares)) and three 
L — 16 (empty squares), 32 (asterisks) and 64 (full 
dimensions, in terms of the reduced variable N = 
We find that P(n) is consistent with a power-law n a 
n (dashed lines). Within numerical accuracy, the 
a = 1.25 both in two and three dimensions. 
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strings of cutting-bonds. Thus the failure of Eq. ( J15| ) for 
percolation clusters simply means that these do not be- 
have as strings of cutting bonds do. In other words, for 
a near 1, it is not correct to approximate a percolation 
cluster as a string of cutting bonds. 



IV. THE ROLE OF BLOBS 

A. Structural fluctuations 

Our results (Fig. |l|) show that the maxflow exponent 
v(a) follows Eq. (Jfj), although fluctuations in the num- 
ber n of cutting bonds, which exist and are relevant in 
real percolation clusters, were disregarded in its deriva- 
tion. So we face a somehow paradoxical situation, since a 



T3, 

I 

> 

1Ef 



I 

J 1.5 

25, 

~ 1 

0.5 



-8 -7 -6 -5 -4 -3 -2-10 1 

a 

FIG. 3: Maxflow (squares) and ohmic current (circles) scal- 
ing exponents, as numerically estimated for strings of cut- 
ting bonds. The number of bonds n on the string is dis- 
tributed according to Eq. (Q), with a — 1.00 and n* L — L ly/ " 
with 1/v = 0.75. Averages were taken over 10 7 samples, for 
L = 32, 128, 512, 2048 and 8192. Notice that (apart from a 
trivial shift) both critical exponents saturate to (a + l)/v for 
large a. For the maxflow exponent this is the behavior pre- 
dicted by Eq. (|l5|). The fact that the (shifted) conductivity 
exponent has the same behavior for all a > indicates that 
the sum of resistances along the string is dominated by the 
largest one, in this regime. 



naive calculation gives the correct result (Eq. (0)), while 
a seemingly more careful calculation that takes into ac- 
count the fluctuations in n (Eq. fll5|)) does not. As men- 
tioned in the previous section, this means that our as- 
sumption that the maximum flow is determined by the 
cutting bonds alone needs to be revised. In order to test 
this assumption, we separately measure the maximum 
flow allowed by cutting bonds and by blobs, which we 
call m c and mj, respectively, for each percolation clus- 
ter. The overall maximum flow is the minimum of these. 
The procedure works such that one picks first the small- 
est of the cutting bond capacities, and then assigns to 
them an infinite capacity. Then the maxflow is found, 
which is now given by the minimal blob min-cut (config- 
uration) . Figures [| a to c show how the PDF of m c 
(cutting-bond flow) and the total PDF vary with a. For 
non-anomalous values a — (Fig. ||a), the distribution 
is centered around a well-defined mean value. With in- 
creasing a one enters the anomalous regime, and the PDF 
develops a power-law tail. This would be expected to re- 
sult from the cutting bonds, while the blob flows mi, have 
a much narrower distribution, decaying roughly exponen- 
tially for large flows. This means that, when m c is large, 
most probably mt will be much smaller and thus the 
overall flow will be determined by nib- Thus, although 
our derivation of Eq. ([l5]) is correct for strings of cutting 
bonds, it is the blobs that determine the flow in those 
rare cases in which m c is large. Therefore the power-law 
tail in P(m c ), which is responsible for the saturation of 
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FIG. 4: Probability distribution for the maxflow allowed by 
cutting bonds (squares) , blobs (circles) and resulting maxflow 
(crosses), which is the minimum of both. Results are shown 
for L = 256 in two dimensions. From top to bottom, the 
disorder exponent is: a — 0.0,0.5 and 0.7. 



v(a) at large values of a in Eq. (|15), is suppressed by 
blobs on percolation clusters. Figure | illustrates this by 
showing that the fraction of cases - for a given maxflow 
m - that are dominated by the blob contribution follow 
a separate PDF. The collapse is not completely perfect, 
since there may be a very slight trend in the total fraction 
of blob-dominated cases with increasing L. On the other 
hand, the variances of the maxflow distributions scale as 
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FIG. 6: The number of SCB's for L = 256 vs. the backbone 
mass, sample-to-sample. 10000 samples. 



expected (as the mean). It is worth mentioning that the 
distribution of fc-cuts (number of bonds in the min-cut) 
is roughly exponential, so that (k) is of the order of 1.4 
... 1.5 for a = 0.7. 

It is also worth pointing out that there are cross- 
correlations between the structural quantities on one 
hand, and between the structure and the maxflow on the 
other hand. These are illustrated in Figures | and g In 
a system with a given L it is after a moment's delibera- 
tion rather clear that there may be an inverse correlation 
between the number of cutting bonds and the samplc- 
to-sample weight of the backbone. We have not tried to 
measure this relationship quantitatively, but given such 
a relation it is no surprise (Fig. fj]) that the mass of the 
backbone correlates strongly with the maxflow value. 
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FIG. 7: The maxflow, for L = 256 and m = -0.7, vs. the 
actual backbone mass. The average is a running average over 
20 samples, with consecutive masses from 10 4 samples. 



B. Blob dominance 

In order to prove that our hypothesis, namely that the 
blobs set the maxflow scale, is correct, we still have to 
show that the blob flow mj, has the right scaling prop- 
erties, i.e. nib ~ L -1 /^ 1 ^"). A complete calculation of 
the maximum flow allowed by blobs would require de- 
tailed information about the blob's internal structure. 
However an estimate can be obtained from the following 
arguments. It is known that the backbone at p c has a hi- 
erarchical, or self-similar, structure jt5). At the top level 
of this hierarchy, the backbone itself can be thought of as 
a string of singly connected (cutting) bonds interspersed 
with blobs. Blobs in turn are loops made of doubly con- 
nected bonds interspersed with smaller blobs and so on, 
as depicted in Fig. || This hierarchical structure has its 
counterpart in a similar classification of surfaces which 
separate the backbone into two pieces (cuts). At the 
top level of this hierarchy are the surfaces {Si } that cut 
the backbone at just one bond, next come those surfaces 
{S 2 } that cut the backbone at exactly two bonds, etc. 
The capacity C(S) of a cut S is defined as the sum of 
the capacities i c of the bonds crossed by it. Because of 
the maxflow-mincut theorem, the maximum flow equals 
the minimum of the cuts' capacities. Our assumption 
that cutting bonds alone determine the maximum flow 
is equivalent to minimizing the capacities among the Si 
alone. We now describe how the next level S2 in this hier- 
archy can be analyzed. Coniglio jL5| has shown that the 
derivative of the spanning probability p' with respect to 
p is proportional to the average number of cutting bonds 
< n >. An extension of his reasoning, due to Kantor [ p0| , 
allows one to write the second derivative of p'(p) with re- 
spect to p at p c as d 2 p' /dp 2 \ Pc ~< n(n — 1) — 2N 2 > \ Pe , 





FIG. 8: Backbone structure. 



where n is the number of cutting bonds and N 2 is 
the number of pairs of doubly connected bonds. Be- 
cause by definition d 2 p' /dp 2 = at p, 



37 1, one finds 
that 2 < N 2 >=< n 2 > - < n > at p c . Since < n k 
L h l v |Q , we conclude that the typical number of pairs 
of doubly connected bonds at p c is N 2 ~ L 2 I V . However 
this alone is not enough to estimate the typical maximum 
flow allowed by doubly connected bonds, for they might 
be grouped into blobs in different ways. Fortunately the 
total number n% of doubly connected bonds at p c can also 
be calculated p^| , and it turns out to be < n 2 >~ L 1 /" . 
This means that the blob statistics is dominated by one 
large ring of roughly L 1 '" bonds and therefore contain- 
ing a number of pairs of doubly connected bonds which 
is of order L 2 I V . Using this information we can now esti- 
mate the maximum flow allowed by blobs at the level of 
doubly connected bonds. This large blob dominates the 
maximum flow since lesser blobs, located somewhere else 
along the backbone, will allow a larger flow. Thus one 
has to find the maximum flow for two parallel strings, 
each containing L x l v cutting bonds. The typical flow al- 
lowed by each string is of order L" 1 / 1 ^ 1 ^"' and therefore 
the typical maximum flow allowed by doubly-connected 
blobs, which is twice this, is of the right order. 

Our reasoning for doubly connected bonds only consid- 
ers typical cases, i.e. fluctuations in the number of doubly 
connected bonds are disregarded. If a particular cluster, 
in addition to having a small number of cutting bonds, 
also has a small number of doubly connected bonds, then 
the next levels in this hierarchy would be relevant. The 
same sort of reasoning can be used at all levels in the 
hierarchy of cuts, but because the algebra becomes too 
complicate for triply connected bonds already, we did 
not test this in detail. However it seems safe to assume 
that min{C(S k )} ~ Z,-V"(i-") for all k > 2 as well. 
Additionally notice that, in order for the mincut to be 
located at triply connected bonds, it is necessary that 
the numbers of singly and doubly connected bonds be 
simultaneously small, an occurrence which arguably has 
a small probability. 

We then see that, in those rare cases in which the num- 
ber of singly connected bonds is small (they allow a large 
flow), blobs take their role thus limiting the flow to a 
value which is typically of order _L~ 1 / zy ( 1 ~ Q ). This then 
shows that the correct value of the exponent v(a) is given 
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by Eq. (|6j). Similar behavior is of course expected for 
other transport properties, e.g. conductivity, in the limit 
of anomalous distributions of bond strengths. This is so 
because in this limit the resistance of the whole cluster 
is dominated by that of the mincut, where conductivities 
are interpreted as critical currents. 

V. CONCLUSIONS 

In this article we have demonstrated that the transport 
problem on percolation clusters still holds surprises. Our 
findings deny the widespread notion that, in the limit 
of anomalous strength distributions, it is the cutting- 
bonds alone that determine the transport properties. We 
show analytically, and confirm numerically that, if blobs 
could be neglected (because of their allowing a larger 
maxflow than cutting bonds) then the overall system's 
behavior would be strongly dependent on the ensemble 
(the cutting-bond PDF tail exponent). This ensemble- 
dependence would come about because the number of 
cutting bonds has a "broad" distribution extending down 
to zero. However the predicted ensemble-dependence is 
not there, as we show numerically on large two and three- 
dimensional systems. Using scaling arguments we then 
demonstrate that it is in fact the blobs that finally de- 
termine the average maxflow. However, we are forced to 
finish with the paradoxical conclusion that though the ex- 
pected mechanism for the maxflow, namely cutting-bond 
dominance, does not work in the anomalous regime (large 
a), the original cutting-bond estimate for the transport 
exponent is nevertheless restored by the limiting effect of 
the blobs. 
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Appendix - an optimal algorithm: We note that 
the augmenting path -method is better here than in 
the general maxflow problem so-called push/relabel pre- 
flow algorithms pi] enjoy the most popularity). This 
is since the structure of the backbone is essentially 
one-dimensional, the number of augmentations remains 
small, of the order of one. To remind the reader, such an 
algorithm consists of flow augmentations, which are re- 
peated till the min-cut is formed (by a surface of blocked 
bonds) and maxflow is reached. For each augmentation 
one needs to establish a path from the "source" to the 
"sink" , which can be done e. g. . by using shortest- 
distance path methods Q. 

The one-dimensional nature means that the backbone 
can be decomposed into strings of subsequent cutting 
bonds Cj and and blobs separating such strings Bi. Thus 
the structure is equivalent to the one-dimensional series 
. . . CiBiCi+i .... In principle one may thus write a more 
efficient algorithm by abandoning the lattice structure, 
and describing the internal geometry of each Bi sepa- 
rately. Thus an optimal version of the algorithm would 
entail the following steps: i): establish the structure 
(Bi, Ci). ii): Find an augmenting path along the chain, 
across all Bi. Hi): Augment flow, that is: find the 
smallest capacity in the Ci, and the smallest capacity in 
all the Bi. This is f±. iv): if f\ equals the minimal 
cutting bond capacity stop, otherwise augment (subtract 
fi from Ci, and the paths inside Bi). v): update the 
paths inside those Bj, only, where a bond was saturated 
by fi (i = 1 to begin with). Goto Hi). We have used, 
instead, an Euclidean background for the maxflow part, 
since the scaling of the matching program is indeed the 
bottleneck. 
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